Functional Analysis of G6PD Variants Associated With Low G6PD Activity in the All of Us Research Program.

Glucose-6-phosphate dehydrogenase (G6PD) protects red blood cells against oxidative damage through regeneration of NADPH. Individuals with G6PD polymorphisms (variants) that produce an impaired G6PD enzyme are usually asymptomatic, but at risk of hemolytic anemia from oxidative stressors, including certain drugs and foods. Prevention of G6PD deficiency-related hemolytic anemia is achievable through G6PD genetic testing or whole-genome sequencing (WGS) to identify affected individuals who should avoid hemolytic triggers. However, accurately predicting the clinical consequence of G6PD variants is limited by over 800 G6PD variants which remain of uncertain significance. There also remains significant variability in which deficiency-causing variants are included in pharmacogenomic testing arrays across institutions: many panels only include c.202G>A, even though dozens of other variants can also cause G6PD deficiency. Here, we seek to improve G6PD genotype interpretation using data available in the All of Us Research Program and using a yeast functional assay. We confirm that G6PD coding variants are the main contributor to decreased G6PD activity, and that 13% of individuals in the All of Us data with deficiency-causing variants would be missed if only the c.202G>A variant were tested for. We expand clinical interpretation for G6PD variants of uncertain significance; reporting that c.595A>G, known as G6PD Dagua or G6PD Açores, and the newly identified variant c.430C>G, reduce activity sufficiently to lead to G6PD deficiency. We also provide evidence that five missense variants of uncertain significance are unlikely to lead to G6PD deficiency, since they were seen in hemi- or homozygous individuals without a reduction in G6PD activity. We also applied the new WHO guidelines and were able to classify two synonymous variants as WHO class C. We anticipate these results will improve the accuracy, and prompt increased use, of G6PD genetic tests through a more complete clinical interpretation of G6PD variants. As the All of Us data increases from 245,000 to 1 million participants, and additional functional assays are carried out, we expect this research to serve as a template to enable complete characterization of G6PD deficiency genotypes. With an increased number of interpreted variants, genetic testing of G6PD will be more informative for preemptively identifying individuals at risk for drug- or food-induced hemolytic anemia.

Analysis of the Combined Effect of rs699 and rs5051 on Angiotensinogen Expression and Hypertension.

Hypertension (HTN) involves genetic variability in the renin-angiotensin system and characterizing this variability will help advance precision antihypertensive treatments. We previously reported that angiotensinogen (AGT) mRNA is endogenously bound by mir-122-5p and that rs699 A>G significantly decreases reporter mRNA in the functional mirSNP assay PASSPORT-seq. The AGT promoter variant rs5051 C>T is in linkage disequilibrium (LD) with rs699 A>G and increases AGT transcription. We hypothesized that the increased AGT by rs5051 C>T counterbalances AGT decrease by rs699 A>G, and when these variants occur independently, would translate to HTN-related phenotypes. The independent effect of each of these variants is understudied due to their LD, therefore, we used in silico, in vitro, in vivo, and retrospective clinical and biobank analyses to assess HTN and AGT expression phenotypes where rs699 A>G occurs independently from rs5051 C>T. In silico, rs699 A>G is predicted to increase mir-122-5p binding strength by 3%. Mir-eCLIP assay results show that rs699 is 40-45 nucleotides from the strongest microRNA binding site in the AGT mRNA. Unexpectedly, rs699 A>G increases AGT mRNA in a plasmid cDNA HepG2 expression model. GTEx and UK Biobank analyses demonstrate that liver AGT expression and HTN phenotypes were not different when rs699 A>G occurs independently from rs5051 C>T, allowing us to reject the original hypothesis. However, both GTEx and our in vitro experiments suggest rs699 A>G confers cell-type specific effects on AGT mRNA abundance. We found that rs5051 C>T and rs699 A>G significantly associate with systolic blood pressure in Black participants in the UK Biobank, demonstrating a 4-fold larger effect than in White participants. Further studies are warranted to determine if the altered antihypertensive response in Black individuals might be due to rs5051 C>T or rs699 A>G. Studies like this will help clinicians move beyond the use of race as a surrogate for genotype.

Clinically important alterations in pharmacogene expression in histologically severe nonalcoholic fatty liver disease.

Polypharmacy is common in patients with nonalcoholic fatty liver disease (NAFLD) and previous reports suggest that NAFLD is associated with altered drug disposition. This study aims to determine if patients with NAFLD are at risk for altered drug response by characterizing changes in hepatic mRNA expression of genes mediating drug disposition (pharmacogenes) across the histological NAFLD severity spectrum. We utilize RNA-seq for 93 liver biopsies with histologically staged NAFLD Activity Score (NAS), fibrosis stage, and steatohepatitis (NASH). We identify 37 significant pharmacogene-NAFLD severity associations including CYP2C19 downregulation. We chose to validate CYP2C19 due to its actionability in drug prescribing. Meta-analysis of 16 independent studies demonstrate that CYP2C19 is significantly downregulated to 46% in NASH, to 58% in high NAS, and to 43% in severe fibrosis. Our data demonstrate the downregulation of CYP2C19 in NAFLD which supports developing personalized medicine approaches for drugs sensitive to metabolism by the CYP2C19 enzyme.

Functional interpretation, cataloging, and analysis of 1,341 glucose-6-phosphate dehydrogenase variants.

Glucose-6-phosphate dehydrogenase (G6PD) deficiency affects over 500 million individuals who can experience anemia in response to oxidative stressors such as certain foods and drugs. Recently, the World Health Organization (WHO) called for revisiting G6PD variant classification as a priority to implement genetic medicine in low- and middle-income countries. Toward this goal, we sought to collect reports of G6PD variants and provide interpretations. We identified 1,341 G6PD variants in population and clinical databases. Using the ACMG standards and guidelines for the interpretation of sequence variants, we provided interpretations for 268 variants, including 186 variants that were not reported or of uncertain significance in ClinVar, bringing the total number of variants with non-conflicting interpretations to 400. For 414 variants with functional or clinical data, we analyzed associations between activity, stability, and current classification systems, including the new 2022 WHO classification. We corroborated known challenges with classification systems, including phenotypic variation, emphasizing the importance of comparing variant effects across individuals and studies. Biobank data made available by All of Us illustrate the benefit of large-scale sequencing and phenotyping by adding additional support connecting variants to G6PD-deficient anemia. By leveraging available data and interpretation guidelines, we created a repository for information on G6PD variants and nearly doubled the number of variants with clinical interpretations. These tools enable better interpretation of G6PD variants for the implementation of genetic medicine.

Quantification of spatial pharmacogene expression heterogeneity in breast tumors.

BACKGROUND: Chemotherapeutic drug concentrations vary across different regions of tumors and this is thought to be involved in development of chemotherapy resistance. Insufficient drug delivery to some regions of the tumor may be due to spatial differences in expression of genes involved in the disposition, transport, and detoxification of drugs (pharmacogenes). Therefore, in this study, we analyzed the spatial expression of 286 pharmacogenes in six breast cancer tissues using the recently developed Visium spatial transcriptomics platform to (1) determine if these pharmacogenes are expressed heterogeneously across tumor tissue and (2) to determine which pharmacogenes have the most spatial expression heterogeneity. METHODS AND RESULTS: The spatial transcriptomics technology sequences the transcriptome of 55 um diameter barcoded sections (spots) across a tissue sample. We analyzed spatial gene expression profiles of four biobank-sourced breast tumor samples in addition to two breast tumor sample datasets from 10× Genomics. We define heterogeneity as the interquartile range of read counts. Collectively, we identified 8887 spots in tumor regions, 3814 in stroma, 44 in lymphocytes, and 116 in normal regions based on pathologist annotation of the tissues. We showed statistically significant differences in expression of pharmacogenes in tumor regions compared to surrounding non-tumor regions. We also observed that the most heterogeneously expressed genes within tumor regions were involved in reactive oxygen species (ROS) handling and detoxification mechanisms. GPX4, GSTP1, MGST3, SOD1, CYP4Z1, CYB5R3, GSTK1, and NAT1 showed the most heterogeneous expression within tumor regions. CONCLUSIONS: The heterogeneous expression of these pharmacogenes may have important implications for cancer therapy due to their ability to impact drug distribution and efficacy throughout the tumor. Our results suggest that chemoresistance caused by expression of GPX4, GSTP1, MGST3, and SOD1 may be intrinsic, not acquired, since the heterogeneity is not specific to chemotherapy-treated samples or cell type. Additionally, we identified candidate chemoresistance pharmacogenes that can be further tested through focused follow-up studies.

Evaluation of US Food and Drug Administration Drug Label Recommendations for Coadministration of Antivirals and Acid-Reducing Agents.

Coadministration with acid-reducing agents (ARAs), including proton pump inhibitors (PPIs), histamine H2 -receptor antagonists (H2 blockers), and antacids has been demonstrated to reduce antiviral exposure and efficacy. Therefore, it is essential that US Food and Drug Administration (FDA) drug labels include recommendations to manage these drug-drug interactions (DDIs). This investigation analyzed information in FDA drug labels to manage DDIs between ARAs and antivirals approved from 1998 to 2019. To ascertain clinical adoption, we assessed whether FDA label recommendations were incorporated into current antiviral clinical practice guidelines. We identified 82 label recommendations for 43 antiviral approvals. Overall, 56.1% of recommendations were deemed clinically actionable, with the most common actionable management strategies being dose adjustment during coadministration (40.2%) and coadministration not recommended (9.8%). The sources informing DDI recommendations were clinical DDI studies (59.8%) and predictions of altered exposure (40.2%). Antivirals with low aqueous solubility were more likely to have label recommendations and were more commonly investigated using clinical DDI studies (P < 0.01). For recommendations informed by clinical DDI studies, changes in drug exposure were associated with actionable label recommendations (P < 0.01). The frequency of exposure changes in clinical DDI studies was similar across antiviral indications, but exposure changes were numerically higher for antacids (71.4%) relative to PPIs (42.9%) and H2 blockers (28.6%). Of DDI pairs identified within drug labels, 76.8% were included in guidelines, and recommended management strategies were concordant in 90.5% of cases. Our findings demonstrate that current regulatory oversight mostly (but not completely) results in actionable label recommendations to manage DDIs for high-risk antivirals.

Mapping the miRNA-mRNA Interactome in Human Hepatocytes and Identification of Functional mirSNPs in Pharmacogenes.

MiRNAs regulate the expression of hepatic genes involved in pharmacokinetics and pharmacodynamics. Genetic variants affecting miRNA binding (mirSNPs) have been associated with altered drug response, but previously used methods to identify miRNA binding sites and functional mirSNPs in pharmacogenes are indirect and limited by low throughput. We utilized the high-throughput chimeric-eCLIP assay to directly map thousands of miRNA-mRNA interactions and define the miRNA binding sites in primary hepatocytes. We then used the high-throughput PASSPORT-seq assay to functionally test 262 potential mirSNPs with coordinates overlapping the identified miRNA binding sites. Using chimeric-eCLIP, we identified a network of 448 miRNAs that collectively target 11,263 unique genes in primary hepatocytes pooled from 100 donors. Our data provide an extensive map of miRNA binding of each gene, including pharmacogenes, expressed in primary hepatocytes. For example, we identified the hsa-mir-27b-DPYD interaction at a previously validated binding site. A second example is our identification of 19 unique miRNAs that bind to CYP2B6 across 20 putative binding sites on the transcript. Using PASSPORT-seq, we then identified 24 mirSNPs that functionally impacted reporter mRNA levels. To our knowledge, this is the most comprehensive identification of miRNA binding sites in pharmacogenes. Combining chimeric-eCLIP with PASSPORT-seq successfully identified functional mirSNPs in pharmacogenes that may affect transcript levels through altered miRNA binding. These results provide additional insights into potential mechanisms contributing to interindividual variability in drug response.

Life-Threatening Docetaxel Toxicity in a Patient With Reduced-Function CYP3A Variants: A Case Report.

Docetaxel therapy occasionally causes severe and life-threatening toxicities. Some docetaxel toxicities are related to exposure, and inter-individual variability in exposure has been described based on genetic variation and drug-drug interactions that impact docetaxel clearance. Cytochrome P450 3A4 (CYP3A4) and CYP3A5 metabolize docetaxel into inactive metabolites, and this is the primary mode of docetaxel clearance. Supporting their role in these toxicities, increased docetaxel toxicities have been found in patients with reduced- or loss-of-function variants in CYP3A4 and CYP3A5. However, since these variants in CYP3A4 are rare, little is known about the safety of docetaxel in patients who are homozygous for the reduced-function CYP3A4 variants. Here we present a case of life-threatening (grade 4) pneumonitis, dyspnea, and neutropenia resulting from a single dose of docetaxel. This patient was (1) homozygous for CYP3A4*22, which causes reduced expression and is associated with increased docetaxel-related adverse events, (2) heterozygous for CYP3A4*3, a rare reduced-function missense variant, and (3) homozygous for CYP3A5*3, a common loss of function splicing defect that has been associated with increased docetaxel exposure and adverse events. The patient also carried functional variants in other genes involved in docetaxel pharmacokinetics that may have increased his risk of toxicity. We identified one additional CYP3A4*22 homozygote that received docetaxel in our research cohort, and present this case of severe hematological toxicity. Furthermore, the one other CYP3A4*22 homozygous patient we identified from the literature died from docetaxel toxicity. This case report provides further evidence for the need to better understand the impact of germline CYP3A variants in severe docetaxel toxicity and supports using caution when treating patients with docetaxel who have genetic variants resulting in CYP3A poor metabolizer phenotypes.

High-throughput screening system for inhibitors of human Heat Shock Factor 2.

Development of novel anti-cancer drug leads that target regulators of protein homeostasis is a formidable task in modern pharmacology. Finding specific inhibitors of human Heat Shock Factor 1 (hHSF1) has proven to be a challenging task, while screening for inhibitors of human Heat Shock Factor 2 (hHSF2) has never been described. We report the development of a novel system based on an in vivo cell growth restoration assay designed to identify specific inhibitors of human HSF2 in a high-throughput format. This system utilizes a humanized yeast strain in which the master regulator of molecular chaperone genes, yeast HSF, has been replaced with hHSF2 with no detrimental effect on cell growth. This replacement preserves the general regulatory patterns of genes encoding major molecular chaperones including Hsp70 and Hsp90. The controlled overexpression of hHSF2 creates a slow-growth phenotype, which is the basis of the growth restoration assay used for high-throughput screening. The phenotype is most robust when cells are cultured at 25 °C, while incubation at temperatures greater than 30 °C leads to compensation of the phenotype. Overexpression of hHSF2 causes overexpression of molecular chaperones which is a likely cause of the slowed growth. Our assay is characterized by two unique advantages. First, screening takes place in physiologically relevant, in vivo conditions. Second, hits in our screen will be of medically relevant potency, as compounds that completely inhibit hHSF2 function will further inhibit cell growth and therefore will not be scored as hits. This caveat biases our screening system for compounds capable of restoring hHSF2 activity to a physiologically normal level without completely inhibiting this essential system.